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ABSTRACT 

We present results from a study of proto-giobular cluster candidates in the interacting 
galaxy system Arp 284 (NGC 7714/5) using data from the Hubble Space Telescope. 
Previous studies of the Antennae and M51 have suggested that the majority of young 
massive star clusters dissolve within 20 Myr due to mass loss. We use the evolutionary 
synthesis code STARBURSt99 to estimate ages and extinctions for approximately 175 
clusters visible with HST. We also use lower-resolution GALEX and ground-based Ha 
data to estimate the ages of the giant H ii regions in which these clusters are found, 
and compare the Spitzer colours of these H ii regions to those of star forming regions 
in other interacting systems. The ages are also used to aid in the interpretation of 
Chandra X-ray data. 

Clusters in the tidal tails of NGC 7714 are generally found to have ages less than 
20 Myr, though observational limits make the significance of this result uncertain. 
Older clusters, though not numerous, have nearly the same spatial distribution within 
the imaged portion of NGC 7714 as young clusters. The cluster population in the 
bridge connecting the two galaxies appears to be older, but the data in this part of 
the system is too limited to draw firm conclusions. The ages of the giant H ii regions in 
NGC 7714 are generally older than those of their constituent clusters, possibly indicat- 
ing that the young clusters we detect are surrounded by their dispersed predecessors. 
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1 INTRODUCTION 

The study of star formation in interacting galaxies has been 
revolutionized by the Hubble Space Telescope. One of the 
early results to emerge from HST was the discovery of a class 
of bright, blue clusters < 300 Myr old in the p eculiar ellipti- 
cal galaxy NGC 1275 l|Holtzman et al.lll992l ). These young 
massive clusters (YMCs) were also found in the prototypical 
merger remnant NGC 7252, with ages ran ging from 34-500 
Myr, consistent with the age of the merger (jWhitmore et al] 
1 19931 ). Many more examples soon followed, in which it was 
found that the properties of these YMCs implied that they 
might b e young globular clusters, though with h igher abun- 
dances (|Conti fc Vaccalll994l : iHunter et ailll994l '). 

In interacting galaxy systems, in which star forming 
regions are abundant, the high resolution of HST has al- 
lowed the study of star cluster demographics. Among the 
results to emerge from such studies is the apparent disso- 
luti on of young star clusters, also known as infant mortal- 
ity. iFall et ahl (|2005l ) used evolutionary synthesis models to 
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estimate the ages of clusters in the Antennae and found 
that 90 per cent of clusters vanish over each age dex from 
lO'' up to 10^ yr, suggesting that cluster di ssolution con- 
tinues over very long time-scales. In contrast. iBastian et al.l 
(2005aj) found in M51 that 70 per cent of clusters dissolve in 
the first 20 Myr, indicating that disruption is concentrated 
on short time-scales. Recent spectroscopic observations of 
star clusters in the Antennae, combined with more careful 
modeling of the star formation history of the system, does 
not show the long-te rm cluster mortality previously claimed 
(IBastian et al.ll2009l ). thus this issue is still controversial. 

Additional evidence for a short dissolution time-scale 
has b een found in comparis ons of cluster and field popula- 
tions. iTremonti et al.l l|200lh obtained UV spectra and evo- 
lutionary synthesis models of both clusters and diffuse re- 
gions of the dwarf starburst galaxy NGC 5253. Unlike the 
clusters, the diffuse region spectra lacked indications of O- 
type stars. It was suggested that stars form entirely within 
clusters, which dissolve on ~ 10 Myr time-scales, preventing 
short-lived O stars from ever entering the field population. 

Early-stage disruption has been ascribed largely to the 
loss of gas left over from star formation, which is swept 



2 B. W. Peterson et al. 



out of the cluster by stellar winds and massive star su- 
pernovae. The resulting mass loss drops the escape veloc- 
ity below the stellar velocities, causing the cluster to ex- 
pand and decrease in surface brightness. The departing 
stars disperse into the field population in ~few 10 ^ yr 
ijBastian fc Goodwinl liooel : iGoodwin fc BastianI l2006l ). On 
longer time-scales (~ 1 0^ yr and longer), interactions with 
giant molecular clouds (jGieles et al.l 12006) and tidal effects 
are expected to become significant disruption mechanisms. 

If YMCs are progenitors of globular clusters, some 
of them must survive their most massive stars, as 
well as external effects th at may lead to dissolution 
(|de Griis fc Parmentieij [2007^ . Interacting galaxy systems 
have high numbers of YMCs, making them ideal targets 
for cluster population studies. To date, too few interacting 
systems have been observed to definitively determine how 
rapidly and by what mechanism massive clusters dissolve. 
Further, both the Antennae and M51 are advanced interac- 
tions with complex star formation histories. It is necessary 
to confirm infant mortality in early-stage interactions where 
the star formation history is simpler, with a particular em- 
phasis on covering a wide range of tidal environments. 

In this paper we focus on the early-stage interacting 
system Arp 284, which at a distance of 38.6 ± 2.7 Mpc 
(1" = 190 pc) according to the NASA/IPAC Extragalac- 
tic Database (NED) has nearly 175 detected star forming 
regions in archival HST images. We use HST colours and 
the evolutionary synthesis code STARBURSt99 (sb99) to ob- 
tain age and reddening estimates for star forming regions. 
We also use lower-resolutio n UV data from th e Galaxy Evo- 
lution Explorer ( G^Z/i^X: iMartin et aLllioOSl ) and ground- 
based Ha to obtain age estimates for larger H ll regions, and 
mid-IR colours from t he Spitzer Space Te l escope Infrared Ar - 
ray Camera (IRAC: IWerner et al.|[20o3 : iFazio et al.ll2004l ) 
to compare the H ll regions to IR emission clumps found 
in other interacting systems. These age results are used to 



aid in the interpretation of the X-ray dat a of ISmith et al 



ll2005al ). which was obtained using Chandra (jWeisskopf et al 
I2OO2I ). 



2 TIDAL STAR FORMING ENVIRONMENTS 
OF ARP 284 

The peculiar morphology of the Arp 284 system results from 
an off-centre collision between two disc galaxies, with the 
primary (NGC 7714) having roughly three times the mass 
of the companion (NGC 7715). The time of closest approach 
was ~ 100 Myr ago. For a discussion of the interaction his- 
tory of the system, th e reader is directed to the models of 
IStruck fc Siiiithl l|2003l ). In this section we summarize previ- 
ous observations of the Arp 284 system that relate to its star 
formation history. Most of the tidal features are labelled in 
Fig. [TJ while the H 11 region complexes in the primary disc 
are labelled in Fig.El 



2.1 Nucleus of NGC 7714 

The NGC 7714 nucleus was the first object to be de- 
scribed using the term 'starburst,' although there are ear- 
lier r eferences in the litera ture to 'bursts of star forma- 
tion' l|Weedman et al.lll98ll ). Strong Ha emission indicates 



a large nuclear Hll region and ongoing star formation 
(|Gonzalez Delgado et alTll995l : ISmith et al.lll997l ). O and B 

stars are required to fit the UV spectrum, with a best fit- 
ting age ~5 Myr and a bur st of star formation favored ov er 
con tinuous star fo r matio n l|Gonzalez Delgado et al.lll999l ). 

iLangon et aD lj200ll ) found that a second, older nuclear 
population is suggested by strong absorption bands of CO in 
the near-IR. In particular, the K-hand spectrum is compat- 
ible with a combination of red supergiants, AGB stars, and 
RGB stars. They concluded that a two-component model is 
probably an oversimplification of the nuclear star formation 
history, preferring a star formation rate that has decreased 
over several 10* yr, either through a smooth exponential de- 
crease or a series of successively smaller bursts. The burst 
model best fits the observed supernova rate with a burst 
15-50 Myr ago, in addition to the most recent burst 5 Myr 
ago. 

Infrared Space Observatory {ISO) and Spitzer mid-IR 
spectra of the nuclear region show strong features from 
polycyclic arom atic hydrocarbons (PAHs) as well as brigh t 
forbidden lines (jO'Halloran et al.ll2000l : iBrandl et al.ll2004l ). 
The metallicity of NGC 7714 peaks in the nucleus, with 
a value near 0.5Z p, and generally decreases w ith distance 
from the nucleus ([Gonzalez Delgado et al.lll995h . 



2.2 Stellar Ring of NGC 7714 (Feature 1) 

NGC 7714 has a partial ring-like structure to the east. This 
feature does not show prominent Ha emission, refiecting 
a lack of ongoing star formation, while the red continuum 
suggests an old stellar population ([Gonzalez Delgado et al.l 
1 19951 : [Smith et al.lll997l ) . Low levels of H I emission indicate 
that there is little hydrogen gas ([Smith et al.lll997l ). 



2.3 Southern Tails of NGC 7714 
(Features 2 and 3) 

Two parallel tidal tails lie to the southwest of the NGC 7714 
nucleus. The outer southern tail, faintly visible in optical 
images, is revealed by 21 cm maps to be part of a large 
Hi l oop which reconnects with the galaxy in the northern 
tail ([Smith fc Wallinlll992l : ISmith et al.lll997l ). The loop is 
also visible in GALEX FUV, which is shown smoothed with 
overlaid H I contours in Fig. |3| This gas loop is an entirely 
separate structure from the inner southern tail, which hosts 
significant ongoing star formation. 

The base of the inner sout hern ta il just south of the nu- 
cleus contains Hll region A (iSmith et al. 2005a). This reg ion 
is vis i ble in 20 cm r a dio co ntinuum ( Gonzalez Delgado et al.l 
[1995I : [Smith et al.l [ 19971 ). The HST images reveal nu- 
merous bright clusters along this structure, with an- 
other large Hll region at the edge of the field (region C; 
[Conzalcz Delgado ct aL 1995), coincident with a bright X- 
ray source ([Smith et al.|[2005ai). The ages of regi o ns A and 
C have been estimated by Garcia- Vargas et al.l ([l997l ) at 
5.0 ± 0.5 and 4.5 ± 0.5 Myr based on photoion ization mod- 
els. T hese regions also appear in Spitzer images ([Smith et al.l 
|2007|V to be discussed below. 

[conzalez Delgado et al.[ ([l995[ ) found Ha region A to 
have a metallicity comparable to the nucleus, about 0.5Zq. 
Region C, at a larger distance from the nucleus, shows a 
lower metallicity of 0.25^q. 
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2.4 Northern Tail of NGC 7714 (Feature 4) 

The northern tail of the primary galaxy is visible but not 
particularly prominent in optical photographs. Near the base 
of the tail, to the northwest of the nucleus, are Hq regions 
B, D, and E. These regions are also detected in 21 cm 
H I emission and red continuum (jSmith et alJ [l993) • HST 
resolves these regions into several smaller groups of clusters. 
The age of region D has been estimated at 3.5 ± 0.5 My r 
using photoionisation modeling arci'a- Vargas et al.ll 19971 ). 

Ha region B has been found to have lower metallicty 
than the nucleus, about 0.25Zq (jConzalez Delgado et al.l 
Il995h . 



2.5 Bridge (Feature 5) 

The bridge consists of a northern and a southern component. 
The northern bridge shows 'clumpy' emission from Hi gas, 
and several Ha peaks are also found, offset slightly to the 
south, indicating ongoing star formation. The southern com- 
ponent lacks H I gas and Ha emission, but is readily seen in 
the red continuum, likely indicating an older stellar popu- 
lation (jSmith et al.lll997r ). No metallicity data is currently 
available for the bridge. 



2.6 Nucleus of NGC 7715 

The lack of significant Ha emission in the NGC 7715 nucleus 
indicates that there is no ongoing massive star formation 
l|Gonzalez Delgado et~ai1ll995l : ISmith et al.lll997l '). Emission 
from the 21 cm Hi radio line is also minimal, and offset to 
the north of the optical galaxy. The nucleus is comparatively 
bright i n the red continuum , suggesting an aged stellar pop- 
ulation (|Smith et al.lll997h . Spectral studies of the nucleus 
show a post-starburst spe ctrum, again co nfirming the lack 
of ongoing star formation l|Bernlohij|l993l ). 



2.7 Eastern Tail of NGC 7715 
(Feature 6) 

The eastern tail of NGC 7715 is faintly seen in optical im- 
ages. Radi o maps show an Hi counterpart, offset s lightly to 
the north (|Smith fc Wallinlll992l : [Smith et al.lll997l '). Massive 
star formation is thought to be weak based on the lack of Ha 
emission. This feature has not been studied in detail, and it 
falls outside the field of the HST images. This poor cover- 
age is unfortunate, because the models suggest that this tail 
curves behind the NGC 7715 disc, passes behind the bridge, 
and comes back into view as the f aint outer southern tail of 
NGC 7714 ('stru ck fc Smithll2003 h . In this scenario, the tail 
contains a signi ficant fraction of the o riginal mass from the 
NGC 7715 disc l|Struck fc Smithll2003D . 



3 OBSERVATIONS 
3.1 HST observations 

The HST Wide Field and Planetary Camera 2 (WFPC2) 
has observed the Arp 284 system through several observ- 
ing programs between 1995 May 13 and 2001 July 10. These 
programs targeted either the starburst nucleus or SN 1999dn 



in the disc of NGC 7714. The disparate nature of the ob- 
servations means that many of the bright clusters are not 
visible in all filters, and in some cases the clusters were not 
found on the same chip across all filters. The field of view is 
36" X 36" on the PC chip and 80" x 80" for the WF chips, 
with pixel sizes of 0.05" pixel"^ on the PC chip and 0.1" 
pixel" ^ on the WF chips. The images and exposure times 
are summarized in Table [1] Note that the number in the 
filter name approximately corresponds to the central wave- 
length of the filter in nanometers. The available filters were 
F380W (~ U), F555W (~ V), F606W (~ R), and F814W 
(~ I)- 

The NGC 7714 nucleus and most of the inner southern 
tail are visible on the PC chip in all filters. One image in both 
F380W and F814W had to be discarded due to excessive 
cosmic ray contamination. This left three good images in 
the F380W filter, two in F555W, one in F606W, and three 
in F814W. In F814W, H ll region C was visible in two images 
on the WF4 chip. 

The entire northern tail was visible on the PC chip in 
the F380W and F606W filters. In F555W and F814W, only 
12 clusters were visible on the PC chip, with the remainder 
visible on WF4 in F814W. 

The bridge was on the seam between the WF3 and WF4 
chips in all four F380W images, so some of the clusters could 
not be measured in this band. In the other filters, it was 
always on the WF3 chip, with two images in F555W, one in 
F606W, and two in F814W. 



3.2 Hll region observations 

Archival GALEX and Spitzer IRAC data are available for 
the Arp 284 system. The GALEX images increase the num- 
ber of colours that can be fit with Sb99 models, while Spitzer 
colours can be compared to those of star forming regions 
in other interacting systems. These instruments have sig- 
nificantly lower resolution than HST, but still allow the 
H II regions A-E to be resolved, as shown in Fig. O Each 
of these regions contains numerous star clusters resolved by 
HST. For clarity, the point-like sources in the HST data will 
be called 'clusters' while these larger structures will be called 
' H II regions.' With Ha luminosities of 10^® — 10*"^ erg s""^ 
(|Smith et al.ll2005al). these a re properly classified as 'giant 
Hll regions' (e.g. lKennicutd[l984l . 

The Hll regions were located in an Ha map, which was 
then registered to the other images. The Ha and accompa- 
nying 7?-band images were obtained using the 1.8 m Perkins 
Telescope of the Ohio State University, using a CCD with a 
pixel size of 0.49" pixel~^. T hese data were published pre- 
viously in ISmith et al.l l| 19971 ). 

GALEX observations are available in both the FUV and 
NUV bands, with effective wavelengths of 1516 and 2267 A. 
The pixel size is 1.5" pixel"'^, which makes even these rela- 
tively large H ll regions difficult to identify. The integration 
times are 4736 s in FUV and 520 s in NUV. 

IR data from Spitzer are available for both the IRAC 
and Multiband Imaging Photometer for Spitzer (MIPS) in- 
struments. The resolution in the MIPS bands is far too low 
for cluster studies, and there are pronounced artifacts due 
to the point spread function because of the brightness of the 
nucleus, so only the broadband IRAC bands (3.6, 4.5, 5.8, 
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and S.Ovim) were used. The data re duction and mosai cins 
for the Spitzer data are described in ISmith et al.l (|2007l ). 



4 HST DATA ANALYSIS 

4.1 Data Reduction 

All of the images were processed by the standard Hubble 
pipeline. The small number of exposures available in some 
data sets made cosmic ray removal a non-trivial task, since 
median combining two images is ineffective. First, we used 
the irae{3 task cosmicrays to remove events a few pixels in 
size that were well above background. Next, images were 
blinked to identify remaining cosmic rays, which were re- 
moved manually using the fixpix task. Finally, images were 
combined using the imcombine task. For F606W, only one 
image was available, so other filters were used for compari- 
son. 

4.2 Cluster Selection 

Clusters were selected using the iraf task daofind. The 
background level a was determined using imstat in a rel- 
atively dark part of the image. The nuclear region is much 
brighter than the rest of the disc, so it was necessary to use 
a higher background for cluster detection, measured at the 
outer edge of the nucleus. We used daofind in each filter 
and then merged the lists to obtain the final sample. Clus- 
ters detected in only two filters were not included in the 
sample. 

We detected 16 clusters in the nucleus at the 5(t level. 
The nucleus is quite crowded, so no attempt was made to 
push the detection threshold any lower. The lowest detected 
flux was 2.28 x 10~" erg s"^ cm"^ A"^ in F814W. 

In the southern tail, 60 clusters were detected at the 5(t 
level and another 34 were detected at the 3(t level. These 
levels correspond to fluxes of 2.17 and 1.32 x 10^^" erg s^^ 
cm~'^ respectively, in F814W. Before using the soft- 

ware, clusters in this part of the system were selected by 
eye. The 3(t set matched well with the results of the visual 
inspection, and were retained for the analysis. 

In the northern tail, daofind detected 60 clusters at the 
5a level, with no additional clusters detected at the 3(t level. 
This is probably due to crowding in the major complexes, 
which resulted in multiple detections on single clusters at 
lower detection thresholds. These clusters may in fact repre- 
sent multiple unresolved star clusters. In a few cases, multi- 
ple resolved clusters had to be measured in a single aperture 
because they were too close together to measure without 
aperture overlap. 

Cluster selection was difficult in the bridge due to the 
higher background level in this part of the system and the 
lower resolution of the WF3 chip compared to the PC chip. 
For these reasons, it was necessary to lower the detection 
threshold in the bridge to 2a, allowing the detection of 34 
clusters. The faintest flux detected in F814W was 2.17 x 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



10~^^ erg s~^ cm~^ A~^. Note that this is higher than the 
3(T level on the PC chip. Only 23 of the clusters could be 
measured in F380W. 

4.3 Photometry of clusters 

Photometry was performed using the iraf task phot in the 
DAOPHOT package, with magnitudes calibrated to the Vega 
magnitude system using the zero point values given in the 
WFPC2 Data Handbook. The magnitudes and ages for the 
cluster sample are given in Table (2] while the extinctions 
and mass estimates are given in Tabled The determination 
of the ages and extinctions will be discussed in § 15.21 The 
mass determination will be discussed in § 16.21 

The tidal tails generally fell on the PC chip. For these 
clusters, we used aperture radii of 3 or 4 pixels for single 
clusters, depending on crowding, and 6 pixels in six cases 
when barely-resolved clusters had to be measured in a sin- 
gle aperture. Sky brightness for background subtraction was 
determined by the mode of the pixel values in an annulus 
of 10-15 pixels, centered on the cluster. In F814W, most of 
the northern tail was on the WF4 chip. These were measured 
using apertures with 2-pixel radii, or 3-pixel apertures for 
the poorly resolved groups, with sky subtraction determined 
by the mode in a 5-8 pixel annulus. These sizes were chosen 
because the WF chips have approximately twice the angular 
size of the PC chip, so the same area of sky was measured 
regardless of which chip a cluster fell on. Very little of the 
northern tail was visible in F555W. 

The nucleus was always on the PC chip, but required 
the use of smaller apertures with radii of 2 or 3 pixels due to 
crowding. Background levels were determined by the mode 
in a 10-15 pixel annulus. 

The bridge was generally on the WF3 chip. In F380W, 
the bridge unfortunately straddled the WF3 and WF4 chips, 
so a few clusters could not be observed in this filter. The 
aperture radius was 3 pixels, with background determined 
by the median in an 8-13 pixel annulus. This aperture size 
was selected because bridge clusters are often isolated, but 
there are several areas where multiple emission sources are 
packed closely together. In such cases, daofind generally de- 
tected only one or two of the brightest sources, which were 
always sufficiently separated to be measured individually. 

Aperture corrections could not be determined on the 
images. It was therefor e necessary to u s e the standard cor- 
rections prescribed bv iHoltzman et al] (Il995l ). The correc- 
tions for WF3 were used for all of the wide field chips. Stan- 
dard corrections are not given for the F380W or F606W fil- 
ters, so corrections for the F336W and F555W, respectively, 
were used. The corrections were typically < 0.3 mag. 

4.4 H II region photometry 

In addition to the five regions A-E, we measured the nucleus 
(designated region F) and four regions (G-I) in the northern 
part of the bridge, shown in Ha and Spitzer 3.6^im in Fig. [l] 
The GALEX and Spitzer images were registered to the HST 
data so that age estimates could be obtained using the total 
flux in the HST bands. This also allowed identification of the 
clusters making up the H ll regions. Most of region C was 
off frame in aU of the HST bands except F814W. Region G 
could not be measured at 4.5ii-m because of a cosmic ray. 
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The aperture radius was selected as 2.45" (5 pixels on 
the Hq image) in all instruments to ensure that the same 
area of sky was represented at all wavelengths. The sky 
annulus widths were 1.82" in HST, 2.45" in Ha, 4.5" in 
GALEX, and 3.6" in Spitzer. The measurements were not 
aperture-corrected because the HST apertures were largely 
empty, and the emission regions were highly pixelated in 
GALEX and Spitzer. 

Table |4] summarizes the photometric measurements. 
The Ha and R b ands were cal i brated to standard units us- 
ing the results of I James et al.l (|2004l ). The ratio of the Ha 
flux to R flux density provided an estimate of the equiva- 
lent width of the Ha line, EW(Ha). This equivalent width is 
included in the Sb99 models and provides another method 
for determining the ages of the emission regions. Ages are 
presented in Table 5 and extinctions in Table [G] These will 
be discussed in § 16.31 



5 HST RESULTS 
5.1 Cluster colours 

Colour-colour diagrams are shown in Fig. [5] In the top left 
panel, we plot the magnitude in the F555W filter minus 
the magnitude in the F606W filter (F555W-F606W) ver- 
sus (F606W— F814W) for all detected clusters for which 
these colours were available. The curves show Sb99 instan- 
taneous star formation models with 0.2Zq metallicity for 
E{B — V) = (solid curve), 0.24 (dashed curve), and 0.5 
(dot-dashed curve) for ages from 10^-10^° yr. The youngest 
ages are at the lowest values of (F606W— F814W). The mod- 
els are shown to give a sense of the effect of reddening, and it 
should be noted that 0.2Zq is not the observed metallicity 
for the nucleus and parts of the southern tail. The bridge 
population is generally bluer in (F555W— F606W) and red- 
der in (F606W-F814W) than the populations of the other 
regions. The southern tail has the opposite tendency, lean- 
ing red in (F555W-F606W) and blue in (F606W-F814W). 
Only a few clusters in the northern tail were on frame in 
F555W; those that were tend toward more intermediate val- 
ues in both colours. The nuclear population is similar. 

The top right panel shows all detected clusters in 
(F380W-F606W) versus (F606W-F814W), along with 
Sb99 models. The youngest ages on the curves appear at 
the lowest values of (F606W-F814W). The middle and bot- 
tom panels show the same colours, with the different regions 
plotted separately and error bars shown. The differences be- 
tween the populations are less striking in this colour combi- 
nation than the one discussed above, and arise primarily in 
(F606W-F814W). 

The median statistical uncertainties in the northern 
tail are 0.081 and 0.095 mag for (F380W-F606W) and 
(F606W— F814W), respectively. In the southern tail, the me- 
dian uncertainties are 0.153 and 0.107 mag; in the nucleus, 
0.050 and 0.040 mag; and in the bridge, 0.167 and 0.123 
mag. The uncertainties are large in the southern tail in part 
because it includes clumps detected at the 3a" level. If only 
5a level clumps are included, the median uncertainties are 
0.114 in (F380W-F606W) and 0.088 in (F606W-F814W). 



5.2 Age determination 

Age and reddening estimates were obtained using the STAR- 
BURST99 v5.1 evolutionary synth esis code l|Leitherer et al.l 
1 19991 : 1 Vazquez fc Leithered [20051 ) . We modeled the colours 
of the clusters ass uming a s ingle, instantaneous burst of star 
formation with a iKroupal (|2002l ) initial mass function over 
the range 0.1-100Mq. Spectral energy distributions (SEDs) 
were calculated for cluster ages of 1 Myr to 10 Gyr and in- 
cluded contributions from Ha emission. The SEDs were red- 
dened from 0- 2 mag of E{B — V) using the reddening law of 
ICalzetti et al.l |l99J)- These model spectra were convolved 
with the bandpasses of the four HST filters. The model age 
and extinction providing the best fit to the observations was 
determined using minimization 

Uncertaint ies in both th e ages and extinctions were de- 
termined as in ISmith et al.l (|2008l ). and are driven largely 
by the degeneracy between age and reddening. Breaking this 
degeneracy is particularly difficult for clusters with coverage 
in only three filters. Such clusters generally have larger age 
uncertainties than those with that were covered in all four 
filters. 

Models were run for three different metallicities: Zq, 
OAZq, and 0.2Zq. The OAZq model is representative of 
the spectroscopically measured abundances in the nucleus 
and Hll region A, while the 0.2Zq model is close to that 
of Hll regions B and C. The solar composition models were 
used primarily for comparison. 

Several colour combinations were used because not all 
clusters were available in all colours. Age estimates obtained 
from Sb99 have been found to be most reliable when (1) 
large numbers of colours are used and (2) the bluest avail- 
able combination of c olours are used. U co verage is partic- 
ularly important (e.g. iHancock et al.ll2008i l. We adhered to 
these results in determining which of the available colour 
combinations would provide the best age estimate for each 
cluster at each metallicity. 

Our best-fitting ages for the clusters are given in Ta- 
ble [21 while E{B — V) and masses are shown in Table O 
The age distributions within the features of the system are 
shown in Fig. [B] Finding charts showing the locations of the 
clusters are given in Figs. [71 1121 with circles indicating ages 
^ 10 Myr, diamonds 11-22 Myr, and squares > 22 Myr. 
The apertures used for photometry of the large H ll regions 
are also shown. The clusters are generally very young, occu- 
pying the two younger age bins (1-10 and 11-22 Myr). Very 
few clusters older than 22 Myr are found in the tidal tails. 

Ages determined by evolutionary synthesis techniques 
generally have large uncertainties. In our sample, the uncer- 
tainties are typically very close to the ages. The median age 
is 14 Myr, with a median negative uncertainty of 5 Myr and 
a positive uncertainty of 10 Myr. 

5.3 Luminosity function 

The luminosities in the F814W band for the clusters Ai'L{u) 
were calculated using a distance of 39 Mpc (NED), using a 
bandwidth of 7.2 x 10^'' Hz. The luminosity distribution, ex- 
cluding the extremely luminous nuclear population, is shown 
in Fig. 1131 and provides one estimate of our sample complete- 
ness. We assume a simple power law form for the luminosity 
function, defined by n{L)dL oc L'^dL. A best-fitting line is 
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determined after the distribution turns over, giving a slope 
of —1.3 with a completeness limit of VT'^Lq. 

The uncertainty in the slope may be estimated as fol- 
lows. The distribution has a flat top, so there are two rea- 
sonable possibilities for the turnover point. Selecting as the 
completeness limit l{f '^LQ, we obtain a slope of slope —1.1. 
Our best estimate for the slope is — 1.3 ± .2. Cluster lumi- 
nosit y functions typically have a — —2 (e.g. lWhitmore et al.l 
1 19991 ). which gives a slope of —1 when binned logarithmi- 
cally. Our luminosity function is consistent with this trend 
to within the uncertainties. 



6 DISCUSSION 

6.1 HST cluster ages 

The m etallicities determined by iGonzalez Delgado et al.l 
l|l995l) aid in the interpretation of our evolutionary synthesis 
results. In the southern tail, Hll region A has a measured 
metallicity ~ O.SZq. The best age estimates should there- 
fore come from our 0.4Zq models, which show this region 
to be dominated by clusters with ages 10 Myr or less. 

The metallicity of H ll region C at the base of the south- 
ern tail is ~ 0.25^0 ([Gonzalez Delgado et al.lll995^ . Com- 
bined with the general metallicity gradient observed in the 
galaxy, we are led to favor the 0.2Zq model in the western 
part of the southern tail, which shows more clusters older 
than 22 Myr than does the OAZq model. Note, however, 
that part of region C was off frame (Fig. Illf) . so we only 
have age estimates for two of its constituent clusters. 

In the northern tail, the dependence of derived age on 
metallicity is less pronounced. Only five clusters change age 
bins between the 0.2Zq and OAZq models, all of which are 
near the nucleus. Most of the clusters are in H ll regions B 
and D. These regions have measured abundances ~ 0.25Zq , 
so our 0.2Zq models should provide the best age estimates. 
The cluster ages are all less than 23 Myr in these mod- 
els, with about half 10 Myr or younger. This is not unex- 
pected, since this p art of the northern tail is detected in Ha 
l|Smith et al.ll 19971 1 and Spitzer 8^m (|Smith et al.ll2007l ') . 

The metallicity of the nucleus has been measured at 
~ 0.5Zq, so the OAZq model is preferred here. Seven of the 
16 clusters have ages ^ 10 Myr, with the remainder between 
14-21 Myr. The gap between these age populations is prob- 
ably not significant, since the rapid evolution of red giants 
after 10 Myr makes exact age determination exceptionally 
difficult and oft en results in gaps between 10-15 Myr (e.g. 
iFafi et al.ll2005h . The younger population fits well w ith the 
UV spectral results of lConzalez Delgado et al] (|l999t ) , while 
the older population fits with one of the bursts expected by 
iLancon et atl (|200ll ). 

Of the 124 clusters detected in the tidal tails of NGC 
7714, those with ages > 22 Myr number only five in the 
OAZq model and 12 in the 0.2Zq model. If an age is as- 
signed to each cluster based on the measured or most likely 
metallicity, there are eight older than 22 Myr. The small 
sample of 16 clusters detected in the nucleus has no clusters 
older than 22 Myr. There is little evidence in NGC 7714 of a 
population older than ~ 20 Myr, despite an interaction that 
had closest approach ~ 100 Myr ago. We also find that the 
cluster age groups are evenly distributed over the galaxy. 



The lack of older clusters detected in the nucleus may be 
due to the difficulty in detecting such clusters amidst the 
high background caused by YMCs. 

These results are similar to those obtained by 
iBastian et"all l|2005bl ') for the interacting galaxy NGC 6872. 
Most of the clusters in NGC 6872 were found by colour fit- 
ting to be younger than 30 Myr. These were spread through- 
out the galaxy, extending far out into the tidal tails. Clusters 
older than 30 Myr were concentrated in the central regions 
of the galaxy. However, the 'central regions' of NGC 6872 
have a spatial size ~ 30 kpc, while the imaged portion of 
the NGC 7714 disc in our data is only ~ 7 kpc. The age dis- 
tribution within NGC 7714 is compatible with that found 
by Bastian et al. in NGC 6872, but a wider field will be 
required to determine if the concentration of old clusters is 
similar in the two galaxies. 

The bridge ages are difficult to interpret because of the 
large colour uncertainties, incomplete coverage in F380W, 
and lack of metallicity data. The metallicity in NGC 7714 
generally decreases with distance from the nucleus. However, 
some of the bridge material may originate in NGC 7715, 
which has unknown metallicity. Nevertheless, a low metal- 
licity is expected due to the abundances measured elsewhere 
in the system and because tidal features in interacting galax- 
ies generally show low abundances. 

Despite the large uncertainties, we note that the only 
evidence for a significant population with age > 100 Myr 
is in the bridge. Out of 34 clusters, the numbers with ages 
> 22 Myr are 12 and 17 in the OAZq and 0.2^0 models, 
respectively. However, only one of these clusters has F380W 
data available in the OAZq case, and nine in the 0.2Zq case. 
The existence of an older population in the bridge is by no 
means certain, but if present could indicate that cluster dis- 
solution mechanisms have acted more weakly here than in 
the tidal tails, which in this system are closer to the nucleus. 
Older clusters should survive more easily in the bridge be- 
cause tidal forces are generally smaller there than in the 
tails. However, the gas expulsion mechanisms thought to be 
responsible for infant mortality operate independently from 
the tidal forces, so the rate at which very young clusters 
dissolve is not expected to change. 

To address the issue of sample completeness, we begin 
with the luminosity function (Fig. I13p . Cutting out clusters 
with luminosities below the completeness limit, we have 49 
left in the tidal tails and only four in the bridge. The bridge 
sample is far too incomplete to permit any reliable analysis. 

Only two clusters in the tidal tails that survive the lumi- 
nosity cut have ages > 22 Myr, and those only in the 0.2Zq 
model. It is unsurprising that older clusters would be prefer- 
entially cut from a luminosity-limited sample, since clusters 
get both redder and fainter as their stellar populations age. 

6.2 HST cluster mass function 

It is common practice in studies of cluster ages to use a mass 
cutoff rather than a luminosity cutoff. Due to evolutionary 
fading, clusters of a given age will remain above the luminos- 
ity limit only if they have a sufficiently high mass. Therefore, 
imposing a luminosity limit will generally under-sample old 
clusters. 

Mass estimates are based on the assumption that clus- 
ters of a given age all have the same mass to light ratio, with 
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luminosity values calculated for clusters with mass 10 Mq. 
The luminosity of a cluster relative to a IO^A/q cluster of 
the same age then gives an estimate of the mass. The un- 
certainties in the masses are quite large, up to a factor ~ 2. 
Mass estimates based on the best-fitting age and extinction 
are shown in Table [3] 

The lowest mass for which the sample is complete was 
determined assuming a simple power law form for the cluster 
mass distribution. For OAZq, we obtain a reasonably good 
fit with a slope ~ —1.4 and a limiting mass of IO^'^Mq. The 
quality of the fit is not high due to a lack of clusters in the 
mass bin centered at 10^"Mq, which may reflect the high 
uncertainties in masses determined in this way. The mass 
estimates for the Q.2Zq model produced a slope of —1.5, 
but with a lower mass limit at about W^-^Mq. The mass 
distribution for the 0.2Zq model is shown in Fig. 1141 

The uncertainties for the mass function fits may be es- 
timated using the same technique as was used for the lumi- 
nosity function. For the OAZq model, the next best choice 
of limiting mass would be 10 '"^ Mq, which gives a slope of 
of -1.0. In the 0.2^0 model, both IQ'^'^ Mq and lO^'^M© 
are reasonable choices, producing slopes of —1.8 and —1.1, 
respectively. Using these alternate choices to estimate the 
uncertainties, the OAZq model has a slope — 1.4±0.4 , with 
smaller values being more likely. For the O.2Z0 model, the 
slope is —1.5 ± 0.4. At the limits of the uncertainties, our 
mass functions agree with the typical power l aw index near 
Q = —2, found for example in the Antennae dZhang fc Falll 
1 19991 ) and Magellanic Clouds (jHunter et a"l]|2003h . Due to 
t he smal l size o f our sample, we are unable to test the claim 
of lGieled (|200^ ) that the cluster initial mass function is trun- 
cated at the high mass end. 

The completeness limit in mass, unlike that in luminos- 
ity, is a function of age due to evolutionary fading. Thus 
the turn-over in Fig. 1141 gives an age-averaged completeness 
limit for mass. For younger systems, the mass completeness 
limit is lower than for older systems. The masses are plotted 
against the ages for the 0.2Zq model in Fig. 1151 with the 
NGC 7714 disc clusters as blue circles and bridge clusters 
as red asterisks. The solid black curve shows the mass of 
a cluster with an F840W luminosity equal to our complete- 
ness limit of 10^'^I/Q, as a function of age. At lower ages, we 
are reaching less massive clusters, and are more deficient at 
older ages. For clusters less than 10 Myr old, we are complete 
to about IO^'^Mq, while for lO'' yr we are complete only to 
IQ^-^Mq. The dashed black line indicates the average com- 
pleteness limit of IQ^'^Mq. The significance of evolutionary 
fading is apparent, as all of the clusters with ages ~ 10* yr 
have masses J5 IO^'^Mq. 

Mass cuts alter the age distributions considerably. In 
the OAZq model, seven of the 27 clusters remaining after 
eliminating those with mass < 10^ °Mq are older than 22 
Myr, with the youngest of these 111 Myr old. These account 
for six of the seven clusters remaining in the bridge (only one 
of which has F380W coverage) and one of six in the northern 
tail, but none of the 14 clusters in the southern tail. 

Of the 63 clusters remaining after eliminating those 
with mass < IO^'^Mq in the 0.2Zq model, 31 are in the 
southern tail, 18 in the northern tail, and 14 in the bridge. 
There are 22 clusters with ages > 22 Myr, 12 of them in 
the bridge, of which five lacked coverage in F380W. In the 
southern tail, seven of 31 clusters are older than 22 Myr 



(five of them detected at 3a level), while the northern tail 
has two of 18 clusters older than 22 Myr. The youngest of 
these clusters has an age of 89 Myr. 

The number of clusters falling into each age bin in the 
0.2Zq model is summarized in Table [7] for all measured 
clusters, clusters above the luminosity limit, and clusters 
above the average mass completeness limit. For comparison, 
we have included columns showing the expected numbers 
of clusters assuming a cluster formation rate of 1 cluster 
Myr~^, and constant cluster formation for the last 10^'^ 
(CCF30) and 10*(CCFioo) yr, only including clusters above 
the mass limits. These are limiting cases, based on Fig. 1141 
The observed number of clusters in the mass limit case fall 
between these two limits. Thus with the available data we 
cannot rule out constant cluster formation with no infant 
mortality. More sensitive data are needed to reduce the com- 
pleteness limits and search for evidence of cluster destruc- 
tion. With the current data set, the number of clusters is 
relatively small and the mass cut-off relatively high, thus 
the results on infant mortality are uncertain. 

6.3 H 11 region age estimates 

Ages for the Hll regions were estimated in two ways. The 
colours obtained from large-aperture photometry in the 
GALEX and HST bands were fit to Sb99 models, using the 
methods described in ij 15.21 In regions E and G-J, the FUV 
and NUV measurements provided only upper limits, so these 
bands were not fit to the models. The range of cluster ages 
within the H 11 regions suggests that the star formation his- 
tories are more complex than the instantaneous starbursts 
used to model the individual clusters. We therefore used 
both instantaneous and continuous star formation models 
to estimate the ages of these regions. The true star forma- 
tion history is likely to fall between these two extremes. 

Ages were also estimated using the Ha equivalent 
widths, which were fit to instantaneous burst Sb99 models. 
These data are shown in Table 5, along with the median age 
of the HST clusters detected in each region. Note that the 
HST sample has no clusters in region E (four were detected 
but were covered in only two filters), only two clusters in 
region C, since it is at the edge of the WFPC2 field of view, 
and in the bridge (regions G-J) only one or two clusters were 
found in each region, so this indicator is only useful in re- 
gions A, B, D, and F. When possible, we also compare with 
published ages determined via photoionisation modeling. 

The ages determined using EW(Ha) are generally con- 
sistent with the median age of clusters found in the region, 
with both showing cluster population ages < 15 Myr. These 
estimates are somewhat higher than those in the literature 
(see Table 5), but are within a factor of ~ 3 in most of the 
large aperture photometry cases. 

The ages determined by fitting Sb99 models to the large 
aperture photometry of the Hll regions are much higher 
than those determined by other methods in regions B, D, 
and F. The area within the apertures occupied by clusters 
is relatively small, so it is possible that faint sources outside 
the clusters are making significant contributions to the total 
luminosity of the Hll regions. If these sources are older and 
redder, the age of the region as a whole could be quite differ- 
ent from that of its bright clusters. These sources could be 
older clusters below our detection threshold or the remnants 
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of clusters that have aheady dissolved and entered the field 
population. 

To quantify the possible influence of faint sources on 
the HST apertures, the fraction of luminosity in the F814W 
band for each of the H ll regions /c originating in the mea- 
sured clusters was determined. The /c values are similar 
throughout the NGC 7714 disc, with ~ 80 per cent of the 
light coming from outside the clusters. Unresolved objects 
could significantly impact the colours, which may explain 
the higher age estimates for regions B, D, and F. 

All of the bridge Hll regions had their lowest age es- 
timates provided by colour fitting, while /c varies. Regions 
G and I had the lowest /c values in the system, while H 
and J have the highest. Visual inspection of regions G and 
I indicates several likely clusters that were not selected by 
daofind. Considering the relatively high noise and low res- 
olution in HST bridge data, it is likely that these results 
are products of the data quality. Deeper imaging will be re- 
quired to learn if the bridge hosts a cluster population that 
differs fi-om that of the NGC 7714 disc. 

Extinction estimates are shown in Table IS] along with 
so me values from the liter a ture. The extinction estimates 
of iGonzalez Delgado et al.l (|l995l ) come from the Balmer 
decrement, while the e x tincti on in the nucleus found by 
iGonzalez Delgado et al.l l|l999l ) was determined using evo- 
lutionary synthesis results for UV continuum flux dis- 
tributions, and was smaller than that found using the 
Balmer decrement by a factor ~ 3. Our extinction esti- 
mate for the nucleus is in go od agreement with that of 
IGonzalez Delgado et al.l (|l999l ). while our other estimates 
are higher than the literature values by a factor ~ 2. This 
is not surprising, since the Balmer decrement tends to give 
high er extinction estimates than the continuum by a factor 
~ 2 jPanelli et al.lll988l : [Storchi-Bergmann et al.lll994h . pos- 
sibly because dust in star forming regions with strong neb- 
ular emission is destroyed or r emoved by ionizing radiation, 
stellar winds, or supernovae (jFanelli et al.|[l98g ). Another 
possibility is that the hot stars responsible for the Balmer 
lines are still associated with the dusty molecular clouds in 
which they formed, while the cooler stars that contribute 
to the continuum hav e drifted away from these regions and 
suffer less extinction l|Calzetti et al.|[l994l ). 

6.4 Spitzer colours 

Figs. [TU] and [T7] show as red circles the Spitzer 
[4.5] -[5.8] vs. [3.6] -[4.5] and [5.8] -[8.0] vs. [4.5] -[5.8] 
colours of the Hii regions. These were calculated from 
the flux densities as described in the IRAC Data Hand- 
book. For comparison, we also show emission 'clumps' 
from several o t her i nteracti ng systems, includin g Arp 
24 fCao & We' '2007'), 82 ('Hancock et al.' |2007| ). 107 
ifSmit h ct al. 2005b) , 285 (Smith ct al. 200S) and NGC 
2207/IC 2163 (|Elmegreen et al.ll2006l ), as green crosses. In 
the Arp 285 system, NGC 2856 tail clump 3 and disc clump 1 
were unusual, and are marked in Fig. [16] as a magenta open 
diamond and cyan limit, respectively. The unusually lumi- 
nous clump 'i' in the eastern tail of NGC 2207 is also shown 
separately, as an open blu e diamond. We also d isplay the 
mean colours of field stars of lWhitnev et al] l|2004h (magenta 
open triangle), MOHI stars (M. Cohen 2005, private commu- 
nication; open blue square), quasars (jHatziminaoglou et al.l 



|2005l : red squares), and diffuse galactic dust (jFlagev et al.l 
[200(1; blue crosses). 

Most of the H ll region Spitzer colours fit well with the 
clumps measured in other systems. This is important, since 
the selection criteria are somewhat different. The Arp 284 
regions were selected based on Ha emission, whereas clumps 
in the other systems were selected based on Spitzer emission. 
While the Ha emission usually had a counterpart in Spitzer 
S.Ovim, in some cases the aperture placement would have 
differed slightly if the 8.0p.m emission had been used for 
selection. 

[4.5] — [5.8] is often used as a gauge of the star forma- 
tion rate, with redder colours indica ting higher rates o f 
mass-normalized star formation (e.g. ISmith et all l2005lj ). 
The Hii regions lie between interstellar matter and stars, 
indicating that, as in other interacting systems, they prob- 
ably have contributions from both. However, none of the 
Arp 284 reg i ons is quite as red as the Arp 285 tail clump. 
ISmith et al.l (|2008l ) suggest that the exceptionally red colour 
of this clump is due to strong emission from interstellar mat- 
ter owing to the very young age of the clump, which was 
determined to be 4 Myr using optical colours. Region J has 
the same estimated age and is only slightly bluer in colour, 
suggesting that it may also have numerous hot young stars. 

The [5.8] — [8.0] colours generally fall close to those of 
interstellar matter, which is not surprising since the 5.8 and 
8.0vim bands are expected to be dominated by interstellar 
dust. 

A few of the regions are unusual. Regions B and C are 
unusually red in [3.6] — [4.5], but slightly blue in the other 
two colours. Region F is redder in [3.6] — [4.5] than clumps 
in other systems with similar [4.5] — [5.8] colour. Region H is 
slightly blue in [5.8] -[8.0]. 

The colours of regions B and C suggest excess emis- 
sion at 4.5jxm, perhaps with a slight deficiency at 5.8 and 
8.0^.m, qualitatively similar to low-metallicity dwarf galax- 
ies. In Figs. [18] and 1191 we compare the Spitzer colours of 
Arp 2 84 H ll regions to those o f the dwarf and spiral galax- 
ies of [Smith fc Hancockl (|2009l ). Regions B and C are not 
as red in [3.6] — [4.5] as some the dwarfs, but are generally 
redder than the spirals. This may indicate that they are de- 
ficient in PAHs, though not extremely so, consistent with 
their moderately low metallicities. 

It is interesting that regions B and C stand out in these 
plots, but the other regions do not. Our colour-based age 
estimate for region B is very similar to that of D, and our 
estimates based on EW(Ha) vary little between regions. 
[Garcia- Vargas et al.l (|l997[ ). however, found regions B and 
C to be extremely young. This would imply more nebular 
continuum and Bra emission, which are two of the factors 
that ap pear to contribute to r eddening [3.6] — [4.5] in the 
dwarfs ([Smith fc Hancock[[2009l ). 

The SED for region B is plotted in Fig. [20] We also show 
as black curves the best-fitting 0.2Zq evolutionary synthe- 
sis models from Sb99 for both continuous (top panel) and 
instantaneous burst (bottom panel) star formation models. 
The curves use the best-fitting ages for each model, which 
are 128j;42^^ Myr for continuous star formation and 88lJ^g 
Myr for the instantaneous star formation model. The models 
include both stellar and nebular emission. The solid green 
lines show the same models without nebular emission. The 
models corresponding to la uncertainty (68 per cent con- 
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fidence) are also shown. All of the curves have been scaled 
to match the observations in F606W. The stellar and stel- 
lar+nebular emission models are nearly indistinguishable, 
indicating that nebular emission is relatively unimportant 
for these inferred ages. The large mid-IR excess could result 
from hot dust continuum emission, PAH lines, or an older 
stellar population, either in clusters that have faded below 
our detection limits or in the field. The excess is diminished 
in the continuous star formation model, probably because it 
includes older stars. 

For comparison, the SED of the nucleus is shown in 
Fig. 1211 The top panel shows curves for a continuous star 
formation model with an age of 42llj^ Myr, while the bot- 
tom panel shows an instantaneous starburst with an age of 
222^g Myr. In this plot we h ave also included th e near-IR 
J, H, and K band fluxes from iLan^on et all (|200l[ ). The nu- 
cleus is known to have a complex star formation history, 
likely including stars older than ~ 100 Myr. The mid-IR ex- 
cess is greater in the instantaneous burst model, likely due 
older stars not included in this model. 



6.5 Comparison with Chandra 

6.5.1 Ultra-luminous X-ray sources (ULXs) 

X-ray observations by ISmith et al.l (|2005al ) detected three 
non-nuclear point sources with Lx > 10^^ erg s~^, the stan- 
dard definition of an ultraluminous X-ray (ULX) source. 
Three other point sources lie just below this limit. The na- 
ture of ULXs is not presently understood. In one possible 
scenario, they could be intermediate mass (100 — IOOOMq) 
black holes (IMBH: [Coibert fc Mushotzkvlll999l ). Alterna- 
tively, they may be stellar mass black holes with beam ed 
iKing et al.ll200ll ) or super-Eddington l|Begelmanll2002l ) X- 
ray emission. IMBHs may preferentially fo rm in massive star 
clusters (e.g. IPorteeies Zwart et al.ll2004al ). It takes at least 
4 Myr for a stellar mass black hole/X-ray binary to form, and 
they are expec ted to be commoii in po pulations with ages up 
to « 100 Myr ijRappaport et al.ll2004l ') In contrast, an IMBH 
may f orm in a YMC in ^ 3 Myr (IPortegies Zwart et al.l 
l2004al ). Evolutionary models of the predicted X-ray emis- 
sion from IMBH-|-stellar companion mass-transfer binaries 
indicate that such ULXs should generally have a ges less than 
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As discussed in ISmith et al.l (|2005al ) , the ULX candi- 



l2004ah . 



dates in NGC 7714/5 exist in a range of environments. A 
luminous X-ray point source is located 1.5" away from the 
nucleus, near cluster nOl, which has an estimated age of 14 
Myr. If this association is correct, the young age is more 
consistent with the IMBH scenario. Another possibly young 
ULX candidate is found within Hll region C, which has an 
age of 11 Myr based on EW (Hq) and a spectroscopically de- 
termined age of 4.5 Myr (iGarcia- Varg as et al. 1997). How- 
ever, this second ULX cannot be associated with a specific 
star cluster, since there are several very close together in the 
area. Finally, one ULX is located in the outer western tail of 
NGC 7714, in an area that is undetected in Ha but is bright 
in GALEX images. Three other ULX candidates are neither 
located in H ll regions nor are they close to optically-selected 
star clusters. Given the variety of environments, from this 
small sample we cannot make any strong conclusions about 
the nature of ULXs as a class. 



6.5.2 Extended X-ray emission associated with star 
forming regions 

In addition to point sources, the Chandra map revealed 
extended X -ray em ission associated with some of the 
Hll regions (|Smith et al, 2005ai)- Fig.[2]shows Chandra 0.3- 
8 keV X-ray contours over a colour scale Hq image of NGC 
7714. Extended X-ray emission is found associated with all 
of the disc Hll regions except D. The bridge regions G-J 
also do not have observed diffuse X-ray emission. 

The possible contributors to the extended X-ray emis- 
sion include multiple unresolved high mass X-ray binaries 
(HMXBs), or hot gas from supernovae and stellar winds im- 
pacting the ambient interstellar medium. Based on Local 
Group studies, the HMXB component is expected to be too 
small to account for this emission in the NGC 7714/5 re- 
gions (|Smith et al.ll2005al ). unless there is an excess above 
Local Group galaxies. For the hot gas component, using 
Sb99 models and a constant X-ray production efficiency 
Lx / Lmech, where Lmech is the mechanical luminosity from 
the supernovae and winds, the ratio of the number of Lyman 
continuum photons NLyC to X-ray luminosity NLyc/Lx 
is e xpected to drop off with age for a star forming region 
( Smith et al.ll20Q5al ). However, the one Hll region complex 
in the NGC 7714 disc that is not detected in extended X- 
ray emission, region D, has an EW(Ha)-based age that is 
somewhat older than some of the other regions (Table [7|, 
in spite of having the highest NLyC / Lx . In contrast, region 
E is bright in the X-rays, with a low NlvC / Lx , but has a 
younger age from its EW(Hq). This suggests that the X- 
ray production efficiency may vary with time and/or from 
region to region, and/or the HMXB contribution may be 
larger than expected in some regions. This result is uncer- 
tain, however, due to the relatively large uncertainties on 
our age estimates and the small number of H ll regions in 
our sample. 



7 CONCLUSIONS 

We have used broadband colours to obtain age estimates for 
174 proto-globular cluster candidates in the Arp 284 system. 
The populations detected in the tidal features of NGC 7714 
are generally quite young (< 20 Myr old). However, when 
the sample is limited by a mass cut, the dominance of young 
clusters becomes unclear in the NGC 7714 disc and vanishes 
in the bridge. Due to the small number of clusters above 
the mass cut, this latter result is highly uncertain. Deeper 
imaging is required to better define the cluster mass function 
and push down the completeness limit, particularly in the 
bridge, which has higher photometric uncertainties and less 
F380W coverage than the rest of the system. 

No substantial differences are found between the popu- 
lations within NGC 7714. The clusters in the nucleus do not 
generally differ in age from those in the tidal structures. It 
therefore appears that the early-stage dissolution processes, 
if present, operate via internal processes, rather than local 
environmental effects, but we cannot rule out the possibility 
of weaker dissolution in the bridge. 

Using larger apertures to study Hll region complexes, 
we find in several cases that the region as a whole appears 
to be older than the clusters detected inside it. This sug- 
gests an older, redder population that is unresolved, possibly 
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the remnants of star clusters that have already dissolved or 
which were not detected as discrete sources because of evo- 
lutionary fading. 

The large aperture results suggest that star formation 
has been ongoing in Hll regions B and D, located in the 
northern tail, for a dynamical time-scale, like that in the 
nucleus. In both cases, the star formation history, whether 
continuous or in discrete bursts, is unknown. We also do 
not know if regions B and D are gravitationally bound, or 
simply regions where successive waves or infall events have 
triggered star formation. Region C is similar, but we have 
insufhcient data to comment on its nature. If in the future 
these regions are found to be bound, they could be inter- 
preted as nearby analogs to the h igh redshift disc clumps 
studied bv lElmegreen et al.| (|2009l '). 

The Arp 284 interaction has been in progress for over 
100 Myr, yet we find little evidence for a significant cluster 
population older than 20 Myr. Thus it is possible that clus- 
ters are dissolving on a relatively short time-scale. However, 
due to completeness issues, with the available data we can- 
not definitively confirm the possibility of infant mortality in 
this system. Further, more sensitive observations are needed 
to search for a deficiency of older clusters. 
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Figure 1. Optical photograph of Arp 284 from I ArpI lll966h . with 
numbers labeling the prominent tidal features. North is up and 
east to the left. The field of view is 4.5' x 2.2'. 
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Figure 4. Hll regions in Ho (top) and Spitzer 3.6|am (bottom). 
The aperture radii have an angular size of 2.45", and are also 
shown on the HST images in Figs. [7 til2l The field of view is 
2.4' X 0.9'. 
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Figure 2. NGC 7714 disc indicating H II regions from 
ISmith et all l|2005al '). The colour scale shows Ho. The contours 
are a Chandra 0.3-8 keV x-ray map. 




Figure 3. Smoothed GALEX FUV (colour) overlaid with VLA 
Hi (contours). Note that the western tail loops back to the main 
galaxy, and this ring-like feature is present in both H I and UV. 
The field of view is 5.1' X 3.5'. 
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Figure 5. In the top panels, red circles represent clusters from 
the bridge, black diamonds the southern tail, blue asterisks the 
northern tail, and black squares the nucleus. The curves show 
Sb99 instantaneous star formation models with 0.2Zq metallicity 
for E(B-V) = (solid curve), 0.24 (dashed curve), and 0.5 (dot- 
dashed curve) . These curves give a sense of the effect of reddening, 
but the metallicity is not representative of the nucleus or parts of 
the southern tail. The middle and bottom four panels show the 
same colours as the upper right panel separated by region, with 
error bars. 
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Figure 6. Histograms of cluster ages from 0.2Zq models in the 
southern tail (top left panel), northern tail (top right panel), and 
bridge (bottom left panel), and from the OAZq model in the 
nucleus (bottom right panel). 
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Figure 7. F606W image of NGC 7714 with northern, central, and 
southern areas indicated. These areas are examined more closely 
in Figs. [8 til2l Several H II regions identified in lower-resolution 
data arc also indicated. Note that region C is partially off frame 
toward the bottom right of the box indicating the southern area. 
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Figure 8. F606W image of the northern tail of NGC 7714, with 
ages indicated for all detected clusters for the 0.2Zq models. 
Circles represent ages 1-10 Myr, diamonds 11-22 Myr, squares 
^ 23 Myr. The metallicities of Hll regions B and D have been 
measured at 0.25^0, while the nucleus has metallicity O.SZq. All 
clusters shown are considered part of the northern tail. 
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Figure 9. F606W image of the NGC 7714 nucleus with ages 
from the OAZq model indicated for all detected clusters. Circles 
represent ages 1-10 Myr, diamonds 11-22 Myr, squares ^ 23 Myr. 
The metallicity has been measured at 0.5Zq. 




Figure 10. F606W image of the central region of NGC 7714 with 
ages from the OAZq model indicated for all detected clusters. 
Circles represent ages 1-10 Myr, diamonds 11-22 Myr, squares 
^ 23 Myr. The nucleus has a measured metallicity of 0.5Zq, so 
this model should give the best age estimates. These clusters are 
considered to be part of the southern tail. 
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Figure 11. F606W image of the southern tail of NGC 7714 with 

ages from the 0.2Zq models indicated for all detected clusters. 
Circles represent ages 1-10 Myr, diamonds 11-22 Myr, squares 
^ 23 Myr. Clusters 10 and 15 are inside Hll region C, which has 
a measured metallicity of O.25Z0 




Figure 12. F555W image with all detected clumps in the bridge 
with ages from the 0.2Zq models. Circles represent ages 1-10 
Myr, diamonds 11-22 Myr, squares ^ 23 Myr. The metallicity of 
the bridge has not been measured, but is likely to be low. 
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Figure 13. Luminosity distribution in F814W for all detected 
clusters outside the nucleus. The luminosity is defined as the lu- 
minosity in the F814W band AuL{u), where Au is the band- 
width of 7.2 X 10^^ Hz. The solid red line is a power law fit to 
the distribution after the turnover with slope —1.3 ± 0.3, giving 
a completeness limit of approximately (5.1 ± 0.2) X IO^Lq. 
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Figure 15. Mass vs. Age for the 0.2Zq model. Blue circles rep- 
resent the NGC 7714 disc, while red asterisks are bridge clusters. 
The dashed black line represents the age-averaged mass complete- 
ness limit of 1O*-*M0 from Fig. [14] The solid black curve shows 
the mass of a cluster with an F840W luminosity equal to our 
completeness limit of IO^'^Lq, as a function of age. 
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Figure 14. Cluster mass distributions for the 0.2Zq model. The 
solid red line is a power law fit to the part of the distribution 
after the turnover. The fit has slope —1.5 ± 0.4 with a limiting 
mass ~ 1O*-*M0. 
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Figure 16. Spitzer [4.5] — [5.8] vs. [3.6] — [4.5] colour-colour plot. 
The Arp 284 Hll regions (open red cir cles) are shown along 
with emission 'clum ps' f rom Arp 24 llCao fc Wd l2007t), 82 
llHancock et al.ll2007f) . 107 llSmith et al."2005b"'). 285 (Sm ith et al.l 
l2008t) . and NGC 2207/IC 2163 (Elmcgrccn ct al. 200^; green 
crosses). Shown separately are Arp 285 NGC 2856 tail clump 
3 (magenta open diamond) and disc clump 1 (cyan limit), as 
well as NGC 2207/IC 2163 'clump i' shown separately (open 
blue diamond) . Also shown are the mean colours of field stars of 
IWhitnev et al.l | |2004) (magenta open triangle), MOIII stars (M. 
Cohen 2005, private communication; open blue square), quasars 
llHat ziminaoglou ct al. 2005; red squares), and diffuse dust in the 
Milky Way lIFlaeev et al.ll2006i ; blue crosses). 
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Figure 18. Spitzer [4.5] — [5.8] vs. [3.6] — [4.5] colour-colour 
plot showing Arp 284 H II regions (o pen red circles). The 
dwarf galaxies of ISmith fc Hancock! | |2009l ) are shown by 
metallicity, with 12 + log(0/H) < 7.9 (red open trian- 
gles), 7.8 < 12 + log(0/H) < 7.9 (green filled circles), and 
12 -I- log(0/H) > 8.2 (blue open squares), along with a sample of 
'normal' spirals (black crosses). 
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Figure 19. Same as Fig.[l8] but for [5.8]-[8.0] vs. [4.5]-[5.8]. 



Figure 17. Same as Fig.[T6] but for [5.8]-[8.0] vs. [4.5]-[5.8]. 
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Figure 20. The FUV-mid-IR SED of region B. The black curves 
are the best-fitting 0.2Zq evolutionary synthesis models includ- 
ing both stellar and nebular emission. The model in the top panel 
is for continuous star formation with an age 
while the bottom panel shows an instantaneous burst with an age 
of 8at\l Myr. The solid green line plotted from 3— 9.5|jm shows 
the models without nebular emission. The short-dashed red and 
dot-dashed blue curves correspond to the youngest and oldest 
models and the associated best-fitting extinctions. The models 
have been scaled to match the observations in F606W. 
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Figure 21. Same as Fig. 1201 but shows region F with OAZq mod- 
els, with ages of 421^^^ Myr for the continuous star formation 
model and 222^g Myr for the instantaneous burst model. 
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Table 1. HST WFPC2 imaging of NGC 7714/5 



Dataset 


Exposure 


Filter 


H II Regions 


Hll Regions 






(sec) 




in PC chip 


in WF chip 


Comment 


U3GQ0201R 


500 


F380W 


nucleus, A, B, part C,D,E 


bridge 


bridge split WF3/4 


U3GQ0202R 


500 


F380W 


nucleus, A, B, part C,D,E 


bridge 


bridge split WF3/4 


U3GQ0203R 


400 


F380W 


nucleus, A, B, part C,D,E 


bridge 


CR on PC chip 












bridge split WF3/4 


U3GQ0204M 


400 


F380W 


nucleus,A,B,part C,D,E 


bridge 


bridge split WF3/4 


U6A01601R 


350 


F555W 


nucleus,A,E 


bridge 




U6A01602R 


350 


F555W 


nucleus, A, E 


bridge 




U2E68801T 


500 


F606W 


nucleus,A,B,part C,D,E 


bridge 




U6A01701R 


350 


F814W 


nucleus, A 


B,C,part D 




U6A01702R 


350 


F814W 


nucleus, A 


B,C,part D 




U6A01801R 


350 


F814W 


nucleus, A, E 


bridge 




U6A01802R 


350 


F814W 


nucleus, A, E 


bridge 


CR on PC chip 



Table 2. Cluster magnitudes and age estimates. The full table is 
available in the electronic edition of the journal. 



ID 


F380W 


F555W 


F606W 


F814W 


Age (Zq) 


Age (0.4Zq) 


Age (0.2Zq) 




mag 


mag 


mag 


mag 


Myr 


Myr 


Myr 


1 


22.46±0.06 




22.64±0.03 


22.18±0.05 




151^6 


15+? 


2 


23.30±0.14 




22.84±0.06 


23.10±0.15 




5l^ 


5lt 


3 


22.33±0.07 




21.23±0.03 


21.67±0.07 


111 


9+2 


3l^ 


4 


22.46±0.08 




22.02it0.04 


22.09±0.08 




611 


7+8 
'-2 


5 


21.25±0.04 




20.86±0.02 


21.06±0.04 




511 


611 


6 


22.18±0.06 




22.04±0.03 


22.28±0.09 


51^ 


511 


611 


7 


21.00±0.03 




21.30±0.02 


21.07±0.03 


511 


14+1 


1511 


8 


22.80±0.11 




21.78±0.04 


21.89±0.05 


311 


411 


41? 


9 


20.84±0.03 




20.82it0.04 


20.74±0.07 


511 


141^ 


1511 


10 


21.13±0.04 




21.43±0.03 


21.07±0.05 


1511 


141^ 


1511 


11 
12 
13 


22.89±0.11 
23.39±0.17 
21.41±0.05 




22.27±0.05 
22.97±0.09 
21.51±0.03 


21.94±0.13 
22.71±0.16 
21.17±0.07 


09+94 
■^"^-30 

611 


— 15 

20lS 
141^ 


9n+102 

20lS 
151^ 


14 


22.30±0.10 




21.62±0.07 


21.89±0.22 


0+3 
•5-2 


4l| 




15 


21.89±0.07 




21.73±0.04 


21.04±0.06 


15lf 


2llt 


20ll 
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Table 3. Cluster E{B — V) and mass estimates. The full table is 
available in the electronic edition of the journal. 



ID E{B - V) 


E{B - V) 


E{B - V) 


M{10'^Mq) 


M(1O6M0) 


M{10'^Mq) 


% 


QAZq 


0.2Zq 


% 


0.4Zq 


0.2Zq 



I 


Q + O 


0+0 


n+0.04 


0.029 


0.029 


0.029 


2 




^^•-^^-0.12 


12+"-'2 


0.012 


0.012 


0.013 


3 




22+"-"' 

^•^^-0.08 


24+"-"6 


0.047 


0.053 


0.061 


4 


n+0. 18 


n+0. 06 

^-0 


n+0. IS 

^-0 


0.013 


0.013 


0.013 


5 


06+° "^ 


10+""'' 

^•J^"_0.04 


06+" "'* 


0.042 


0.074 


0.065 


6 


f,+0.08 


n+0.04 

^-0 


n+0.04 

"-0 


0.011 


0.017 


0.017 


7 


0+0 
"-0 


0+0 
^-0 


0+0 

^-0 


0.081 


0.074 


0.051 


8 




n qo + 0.08 


34+0.08 


0.075 


0.086 


0.067 


9 


0+0 


0+0 

"-0 


0+0 

"-0 


0.110 


0.100 


0.070 


10 


0+0 


0+0 


0+° 


0.081 


0.074 


0.081 


11 


n+0. 24 
^-0 






0.063 


0.067 


0.074 


12 


n+0. 22 


0+0.22 
"-0 


n+0. 26 

"-0 


0.024 


0.024 


0.039 


13 


0+0 


0+0 
"-0 


0+0 
"-0 


0.074 


0.067 


0.031 


14 




n 10+0.16 


22+"-l2 


0.053 


0.043 


0.031 


15 


n+0. 06 


06+° "* 


08+° "** 


0.145 


0.142 


0.083 



22 B. W. Peterson et al. 



Table 4. H II region photometry 



Region 


FUV^ 


NUV" 


F380W'' 


F555W" 


F606W'' 


F814W'' 




3.6 nm^ 






A 


2080±95 


1890±81 


549 ±4 


396 ± 1 


459 ± 1 


276 ± 1 


139 ± 14 


2.56 ±0.07 


1.58 ±0.04 




B 


1130±110 


474±63 


468 ±4 




202 ±2 


110 ±2 


92.9 ±9.4 


0.431 ±0.010 


0.391 ± 0.006 




C 


1550±80 


551±40 








153 ± 1 


52.8 ±5.3 


0.235 ±0.019 


0.235 ± 0.010 


1 


D 


500±91 


236±58 


258 ±3 




140 ± 1 


86.2 ±1.7 


32.0 ±3.5 


0.357 ±0.061 


0.255 ± 0.048 


1 


E 


< 218±25 


^ 64.8±21.2 


15.8 ±3.0 




11.2 ±0.7 




7.32 ± 0.79 


^ 0.141 ±0.011 


^ 0.0851 ±0.0150 


^ 


F 


13600±100 


6220±110 


13200 ± 100 


9080 ± 10 


8610 ± 10 


5690 ± 10 


1140 ± 110 


14.5 ±0.1 


12.9 ± 0.1 




G 


^ 70.3±11.2 


^ 22.4±8.3 




17.2 ± 0.4 


16.9 ±0.4 


8.32 ±0.28 


2.8 ±0.3 


0.0152 ±0.0034 




0. 


H 


< 93.2±10.4 


< 43.2±6.5 




27.6 ±0.4 


25.2 ±0.4 


10.9 ±0.3 


4.87 ±0.51 


0.0305 ± 0.0035 


0.0193 ± 0.003 


0. 


I 


^ 39.2±7.3 


^ 32.0±5.9 




14.9 ± 0.4 


14.8 ±0.4 


9.73 ±0.27 


2.02 ± 0.23 


0.0259 ± 0.0032 


0.0186 ± 0.003 


0. 


J 


41.0±6.3 


34.8±5.7 




10.4 ±0.4 


13.0 ±0.4 


5.07 ±0.28 


3.50 ±0.36 


0.0211 ±0.0033 


0.0111 ±0.005 


0. 



FUV, NUV, F380W, F555W, F606W, F814W in units of W'^^ erg s'^ cra'^ A" 
^ Ha in units of 10~^^ erg s~^ cm"^ 
^ 3.6, 4.5, 5.8, 8.0iam in units of mjy 



Table 5. Hll region ages (Myr) 



Region Age (inst^ 


Age (inst)° 


Age (cont)" 


Age (cont)" 


Med. Age_^ Med Age' 


EW(Ha) 


EW(Ha) Age 


EW{Ha) Age Lit. Age_^ 


0.4% 


O.2Z0 


O.4Z0 


O.2Z0 


O.4Z0 ~ O.2Z0 


A 


O.4Z0 


O.2Z0 



A 
B 




10+1 


''oa_207 


128^1272 


5tl 
14tJ 


5tl 
15t^ 


29.8 ±3.3 
39.6 ± 5 


10+; 

g + O 


lltS 

lots 


5 ±0.5 
3.5 ±0.5 


C 














27.3 ±3 


10+; 


11+1 

11-0 
12+i 

9tl 


4.5 ±0.5 


D 
E 








i2mt\%. 


16t^ 


17ti 


22.4 ±3.4 
56.2 ±21 


11+1 

11-0 

8+? 




F 
G 
H 
I 
J 


222+1 
6+2 

y+376 


86+i 
6+2 

O+1092 
°-2 

4ti 


421+^2 

1+15 
^-1 

It? 
1 +3799 


2600tl°g 

1+19 

^-1 

9+8 

^-2 
0+4298 
^-2 

It? 


iota' 

Qi +959 
^^-15 
182+^^7 

13+12 


15tl 
266t^f, 

68t^^ 
12011^4 
5tl 


15 ± 1.6 
28.5 ±7 
27.3 ±4.4 
12.1 ± 2.1 
113 ±82 


12+° 
10+i 

10+; 

13+1 

g+3 


14+; 

iiti 

iitS 
16+; 

7+3 
' -1 


5; 15 - 5( 



Ages in Myr from large-aperture broadband photometry for instantaneous starburst (inst) and continuous star formation (cont) models 

Median age of JJST-resolved clusters in the region. These ar e useful only for regions A. B. D, F; see text for explanation. 

Ag es deter mined by p hotoionization modeling. A, B, C from [Garcia- Vargas et al] |l993); F from [Gonzalez Deleado et al.l (Il995l) and 
iGonza lcz Dcl gado et al I (l999) 

Fract ion of H II region lum inosity contained in sample clusters in F814W 
e From lSmith et al.l ([2005a[') 



Star Clusters in Arp 284 



Table 6. Hll region extinctions 



Region 


E{B - V) (inst) 


E{B - V) (inst) 


E{B - V) (cont) 


E{B - V) (cont) 


Lit. E{B - V)^ 








U.4Z'Q 






A 
B 


r, 99+0.02 
^^■^^-0.02 
n+0.06 


14+002 
^•1^-0.02 

12+0-08 
'-'■1^-0.08 


'^•'^°-0.02 

12+0-08 
^•1^-0.10 


34+0.02 
'-'•'-'^-0.02 

n 10+0. 08 


0.25 
0.28 


C 










0.34 


D 

E 


10+°" 


^-^^^-0.20 


n 1S+0-12 


n+0.48 




F 
G 
H 
I 
J 


p ,0+0.02 

■^■^^-0.24 

"■^^-0.06 
24+°-30 

42+0 12 
'-'■^^-0.12 


r, 90+0.02 

U.28_p 02 

o.28«:;^ 

0.38«:- 


26+"-°^ 

'^•^°-0.02 
22+°-" 

o-44l»:« 
o.i4i«:f2 


0+0 

44+0-44 
^^■^^-0.44 

12+°l^ 
'-'-1^-0.12 


0.22; 0.08 



° A,B , C from lGonzalez Delgado et al.l lll995h : F from I Gonzalez Deleado et al.l lll995h (Balmer decrement) and lGonzalez Delgado et al.l 
lll999l'l (UV continuum) 



Table 7. Clusters by age bin for Z = 0.2Zq 



Age 


All 


Luminosity Cut 


Mass Cut 


CCFg,;! 


CCFicii 


1-10 Myr 


59 


23 


13 


10 


10 


11-22 Myr 


70 


28 


28 


12 


12 


23-100 Myr 


2 





1 


8 


78 


> 100 Myr 


27 


2 


21 








Total 


158 


53 


63 


30 


100 



Expected numbers for constant cluster formation over the last 30 Myr at a rate of 1 cluster Myr 1 
^ As above for 100 Myr 



